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Abstract 

Tensile creep tests, tensile relaxation tests and a tensile test with a constant 
rate of strain are performed on injection- molded isotactic polypropylene at room 
temperature in the vicinity of the yield point. A constitutive model is derived for 
the time-dependent behavior of semi-crystalline polymers. A polymer is treated as 
an equivalent network of chains bridged by permanent junctions. The network is 
modelled as an ensemble of passive meso-regions (with affine nodes) and active meso- 
domains (where junctions slip with respect to their positions in the bulk medium 
with various rates) . The distribution of activation energies for sliding in active meso- 
regions is described by a random energy model. Adjustable parameters in the stress- 
strain relations are found by fitting experimental data. It is demonstrated that the 
concentration of active meso-domains monotonically grows with strain, whereas 
the average potential energy for sliding of junctions and the standard deviation of 
activation energies suffer substantial drops at the yield point. With reference to 
the concept of dual population of crystalline lamellae, these changes in material 
parameters are attributed to transition from breakage of subsidiary (thin) lamellae 
in the sub-yield region to fragmentation of primary (thick) lamellae in the post-yield 
region of deformation. 
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1 Introduction 



This paper is concerned with the experimental study and modelUng of the time-dependent 
behavior of isotactic polypropylene (iPP) at strains up to 20 % in uniaxial tensile tests, 
creep tests and relaxation tests at room temperature. Isotactic polypropylene is chosen 
for the investigation because of its numerous applications in industry (oriented films for 
packaging, reinforcing fibres, nonwoven fabrics, blends with thermoplastic elastomers, 
etc.). 

The nonlinear viscoelastic response of polypropylene was studied by Ward and Wolfe 
(1966) and Smart and Williams (1972) three decades ago, and, in recent years, by Wort- 
mann and Schulz (1994, 1995), Ariyama (1996), Ariyama et al. (1997), Dutta and Edward 
(1997), Read and Tomlins (1997) Tomhns and Read (1998), and Sweeney et al. (1999). 

Viscoplasticity and yielding of iPP have been investigated in the past five years by 
Kalay and Bcvis (1997), Coulon et al. (1998), Seguela et al. (1999), Staniek et al. (1999), 
Nitta and Takayanagi (1999, 2000) and Labour et al. (2001), to mention a few. 

Dynamic mechanical analysis reveals that the loss tangent of iPP demonstrates two 
pronounced maxima being plotted versus temperature (Andreassen, 1999; Seguela et al., 
1999; Lopez- Manchado and Arroyo, 2000). The first maximum (/3-transition in the inter- 
val between T = —20 and T = 10 °C) is associated with the glass transition in a mobile 
part of the amorphous phase, whereas the other maximum (a-transition in the interval 
between T = 70 and T = 110 °C) is attributed to the glass transition in the remaining 
part of the amorphous phase, the so-called "rigid amorphous fraction" (Verma et al., 
1996). This conclusion is confirmed by DSC (differential scanning calorimetry) traces of 
quenched polypropylene which show (in the heating mode) an endotherm near T = 70 °C 
ascribed to thermal activation of amorphous domains with restricted mobility (Seguela et 
al., 1999). 

Isotactic polypropylene is a semi-crystalline polymer containing three different crys- 
tallographic forms: monoclinic a crystallites, hexagonal /? structures, orthorhombic 7 
polymorphs, together with "smectic" mesophase (lijima and Strobl, 2000). At rapid cool- 
ing of the melt (at the stage of injection molding), a crystallites and smectic mesophase 
are mainly developed, whereas /3 and 7 polymorphs are observed as minority components 
(Kalay and Bevis, 1997; Al-Raheil et al., 1998) that disappear after annealing above 130 
°C (Al-Raheil et al., 1998; Labour et al., 2001). 

A unique feature of a spherulites in iPP is the lamellar cross-hatching: development of 
transverse lamellae oriented in the direction perpendicular to the direction of radial lamel- 
lae (lijima and Strobl, 2000; Maiti et al., 2000). The characteristic size of a spherulites 
in injection-molded specimens is estimated as 100 to 200 /xm (Kalay and Bevis, 1997; 
Coulon et al., 1998). These spheruhtes consist of crystaUine lamellae with thickness of 
10 to 20 nm (Coulon et al., 1998; Maiti et al., 2000). The amorphous phase is located 
(i) between spherulites, (ii) inside spherulites, in "liquid pockets" (Verma et al., 1996) 
between lamellar stacks, and (iii) between lamellae in lamellar stacks. It consists of (i) 
relatively mobile chains between spherulites, in liquid pockets and between radial lamellae 
inside lamellar stacks, and (ii) severely restricted chains in the regions bounded by radial 
and tangential lamellae. 

Stretching of iPP specimens results in inter-lamellar separation, rotation and twist of 



2 



lamellae, fine and coarse slip of lamellar blocks and their fragmentation (Aboulfaraj et 
al., 1995; Seguela et al., 1999), chain slip through the crystals, shding and breakage of tie 
chains (Nitta and Takayanagi, 1999, 2000), and activation of rigid amorphous fraction. At 
large strains, these morphological transformations lead to cavitation, formation of fibrills 
and stress-induced crystallization (Zhang et al. 1999). 

It is hard to believe that these mechanically-induced changes in the micro-structure 
of iPP can be adequately described by a constitutive model with a small number of 
adjustable parameters. To develop stress-strain relations, we apply a method of "homog- 
enization of micro-structure," according to which a sophisticated morphology of isotactic 
polypropylene is modelled by an equivalent phase whose deformation captures essential 
features of the response of this semi-crystalline polymer. We choose a network of chains 
as the equivalent phase for the following reasons: 

1. The viscoelastic response of isotactic polypropylene is conventionally associated with 
rearrangement of chains in amorphous regions (Coulon et al., 1998). 

2. Sliding of tie chains along and their detachment from lamellae play the key role in 
the time-dependent response of iPP (Nitta and Takayanagi, 1999, 2000). 

3. The viscoplastic flow in semi-crystalline polymers is assumed to be "initiated in the 
amorphous phase before transitioning into the crystalline phase" (Meyer and Pruitt, 
2001). 

4. The time-dependent behavior of polypropylene is conventionally modelled within the 
concept of a network of macromolecules (Sweeney and Ward, 1995, 1996; Sweeney 
et al., 1999). 

Isotactic polypropylene at room temperature (i.e., above the glass transition tempera- 
ture for the mobile amorphous phase) is treated a permanent network of macromolecules 
bridged by junctions (physical cross-links, entanglements and lamellar blocks). The net- 
work is assumed to be highly heterogeneous (this inhomogeneity is attributed to inter- 
actions between lamellae and surrounding amorphous regions, as well as to local density 
fluctuations in the amorphous phase), and it is thought of as an ensemble of meso- regions 
(MR) with different potential energies. Two types of MRs are distinguished: (i) active 
domains, where junctions can slide with respect to their positions in the bulk material 
as they are thermally agitated (the mobile part of the amorphous phase), and (ii) pas- 
sive domains, where sliding of junctions is prevented by surrounding lamellae (the rigid 
amorphous fraction). 

Straining of a specimen induces 

1. Sliding of meso-domains with respect to each other (which reflects fragmentation 
and coarse slip of lamellae) . 

2. Sliding of junctions with respect to their positions in the stress- free medium (which 
is associated with shp of tie molecules along lamellae and fine slip of lamellar blocks). 

Shding of MRs with respect to each other is modelled as a rate-independent process 
that describes the elastoplastic response of iPP. Sliding of junctions in meso-domains is 
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treated as a rate-dependent phenomenon that reflects the viscoplastic behavior of isotactic 
polypropylene. 

Stretching of a specimen results in an increase in the concentration of active MRs and 
changes in the distribution of meso-domains with various activation energies for sliding 
of junctions driven by release of part of the rigid amorphous fraction due to lamellar 
fragmentation. 

The objective of this study is two-fold: 

1. To report experimental data in a tensile test with a constant strain rate, in creep 
tests and in relaxation tests at several elongation ratios on injection-molded iPP 
specimens annealed for 24 h at the temperature T — 140 °C. 

2. To derive constitutive equations for the time-dependent response of a semi-crystalli- 
ne polymer and to find adjustable parameters in the stress-strain relations by fitting 
observations. 

In previous studies on modelling the response of amorphous and semi-crystalline poly- 
mers in the sub-yield and post-yield regions, see, e.g., pioneering works by Haward and 
Thackray (1968) and G'Sell and Jonas (1979), and more recent publications by Boyce 
et al. (1988), Bordonaro and Krempl (1992), Arruda et al. (1995), Hasan and Boyce 
(1995), Krempl and Bordonaro (1995, 1998), Zhang and Moore (1997), Spathis and Kon- 
tou (1998), Drozdov (2001), Duan et al. (2001), stress-strain curves, creep curves and 
relaxation curves were treated independently of each other (in the sence that different 
adjustable parameters were determined by matching observations in different tests). The 
aim of the present paper is to approximate experimental data in three conventional types 
of mechanical tests within one constitutive model. This allows two approaches in the non- 
linear viscoelasticity with an "internal time" to be compared, the so-called models with 
stress- and strain- induced material clocks (Lustig et al., 1996; Krempl and Bordonaro, 
1998, Wineman, 2002), as well as to shed some light on a mechanism for mechanically- 
induced changes in relaxation (retardation) spectra in the vicinity of the yield point. 

The exposition is organized as follows. The specimens and the experimental procedure 
are described in Section 2. The distribution of meso-regions with various potential energies 
for shding of junctions is introduced in Section 3. Kinetic equations for sliding of MRs 
with respect to each other are proposed in Section 4. Sliding of junctions in active meso- 
domains is modelled in Section 5. The strain energy density of a semi-crystalline polymer 
is determined in Section 6. Constitutive equations arc derived in Section 7 by using the 
laws of thermodynamics. These equations are further simplified to describe observations in 
"rapid" tensile tests, creep tests and relaxation tests. Section 8 is concerned with fitting 
experimental data. Our findings are briefly discussed in Section 9. Some concluding 
remarks are formulated in Section 10. 

2 Experimental procedure 

Isotactic polypropylene (Novolen llOOL) was supplied by BASF (Targor). ASTM dumb- 
bell specimens were injection molded with length 148 mm, width 10 mm and height 3.8 
mm. To erase the influence of thermal history, the samples were annealed in an oven at 
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140 °C for 24 h and slowly cooled by air. To minimize the effect of physical aging on 
the time-dependent response of iPP, mechanical tests were carried out a week after the 
thermal pre-treatment. 

Uniaxial tensile relaxation tests were performed at room temperature on a testing 
machine Instron-5568 equipped with electro-mechanical sensors for the control of longi- 
tudinal strains in the active zone of samples (with the distance 50 mm between clips). 
The tensile force was measured by a standard load cell. The engineering stress, cr, was 
determined as the ratio of the axial force to the cross-sectional area of the specimens in 
the stress- free state. The specimens were loaded with the cross-head speed 5 mm/min 
(that corresponded to the Hencky strain rate tn = 1-1 • 10~^ which provides nearly 
isothermal test conditions (Arruda et al., 1995). 

The engineering stress, cr, is plotted versus the elongation ratio A in Figure 1. The 
true longitudinal stress, Ut, is calculated as Ut = crA (this formula is based on the in- 
compressibility condition). The true stress is also depicted in Figure 1, which shows that 
necking of samples does not occur at elongations up to 50 %. The apparent yield strain, 
eyi, calculated as the intersection point of tangent lines to the true stress-elongation ratio 
diagram at "small" and "large" deformations, equals 0.04. The yield strain, ey2, deter- 
mined as the strain corresponding to the maximum on the engineering stress-engineering 
strain curve, equals 0.08. 

A scries of 8 creep tests was performed at the engineering stresses a? = 10.0 MPa, 

= 15.0 MPa, = 20.0 MPa, = 25.0 MPa, = 30.0 MPa, = 30.38 MPa, 
(T7 = 30.94 MPa, and Ug = 32.80 MPa. The last three values of stress correspond to the 
initial strains (at the beginning of the creep tests) 6° = 0.45, 6° = 0.50, and 6° = 0.60. 
The first four tests were carried out in the sub-yield region of deformation (the initial 
strains are less than the yield strain eyi), whereas the other tests were performed in the 
interval between the yield strains eyi and ey2. 

Each creep test was carried out on a new sample. In the mth test (m = 1, . . . , 8), a 
specimen was loaded with the cross- head speed 5 mm/min up to the engineering stress 
that was preserved constant during the creep test, tc = 20 min. The longitudinal strains, 
e, measured in the first 6 tests are plotted versus the logarithm (log = log^g) of time t 
(the instant t = corresponds to the beginning of a creep test) in Figure 2. This figure 
demonstrates that the rate of increase in strain, e, with time, t, is relatively low at small 
stresses, and it noticeably grows with a. The creep curves at stresses exceeding 20 
MPa are plotted in Figure 3, where the Hencky strain en = In A is presented as a function 
of time t. Figure 3 reveals the primary creep of iPP at the stress a^, transition from the 
primary creep to the secondary creep at the stress cr^, and transition from the secondary 
creep to the ternary creep at higher stresses, cxg to cxg. These transitions are indicated by 
fines AA' (en = 0.06) and BB' (en = 0.14) in Figure 3. 

A scries of 4 relaxation tests was performed at the strains e? = 0.05, = 0.10, 
€3 = 0.15, and = 0.20. The first test was carried out at the strain belonging to the 
interval between the yield strains eyi and ey2, whereas the other tests were performed at 
strains in the post-yield region of deformation. 

Any relaxation test was carried out on a new sample. In the mth test {m — 1, . . . , 4), 
a specimen was loaded with the cross-head speed 5 mm/min up to the longitudinal strain 

that was preserved constant during the relaxation time — 20 min. The engineering 
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stress, cr, is plotted versus the logarithm of time t (the instant t = corresponds to the 
beginning of a relaxation test) in Figure 4. This figure shows that the relaxation curves 
are strongly affected by strain. 



3 Distribution of meso-regions 

A semi-crystalline polymer is treated as a permanent network of chains bridged by junc- 
tions. The network is thought of as an ensemble of meso-regions with various potential 
energies for slippage of junctions with respect to their positions in the reference state. 
Two types of meso-domains are distinguished: passive and active. In passive MRs, all 
nodes arc assumed to move affinelly with the bulk material. In active MRs, the junctions 
slide with respect to their reference positions under loading. 

Sliding of junctions in active MRs is modelled as a thermally activated process. The 
rate of sliding in a MR with potential energy lD is given by the Eyring equation (Eyring, 
1936) 

where is Boltzmann's constant, T is the absolute temperature, and the pre-factor Tq is 
independent of energy lj and temperature T. Confining ourselves to isothermal processes 
and introducing the dimensionless activation energy uj — uj/ik-oTo), where Tq is some 
reference temperature, we arrive at the formula 

r(a;) = roexp(-a;). (1) 

Some lamellae (that restrict mobihty of junctions in passive MRs) break under straining, 

which implies an increase in the concentration of active meso-domains. As a result, the 
number of strands in active MRs grows, and the number of strands in passive meso- 
domains decreases. Denote by A^a(^; ^) the number of strands (per unit mass) in active 
meso-domains with energy uj at instant t > 0. The total number of strands in active MRs, 
Xa(t), reads 

too 

X,{t) = / N,{t,u;)dw. (2) 







The number of strands (per unit mass) in passive MRs, Xp{t), is connected with the 
number of strands in active meso-domains, Xg,(t), by the conservation law 

X,{t)+Xp{t)^X, (3) 

where X is the number of active strands per unit mass (this quantity is assumed to be 
time-independent) . 

The distribution of strands in active MRs is described by the ratio, p{t,uj), of the 
number, Ni^{t,uj), of strands in active meso-domains with energy uu to the total number 
of strands in active MRs, 

iVa(t,Cu) ... 
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and by the concentration, K,{t), of active MRs, 



«(() = (5) 

In what follows, constitutive equations for a semi- crystalline polymer will be derived for 
an arbitrary distribution of active MRs. To fit experimental data, the random energy 
model is employed with 



p{t,uj) =po{t)exp 



^(t))^] (^>o), p{t,uj) = (a;<0). (6) 



2E2(i) 



Here ft is the average activation energy in an ensemble of active meso-domains, E is the 
standard deviation of potential energies for sliding of junctions, and po is determined by 
the condition 

/ p(t,u;)d(jj = 1. (7) 
Jo 

4 Sliding of meso-regions 

It is assumed that meso-domains are not rigidly connected, but can slide with respect 
to each other under straining. Sliding of meso-domains is treated as a rate-independent 
process and is associated with the elastoplastic behavior of a semi-crystalline polymer. 

We suppose that an increase in strain, e, by an increment, de, causes growth of the 
elastoplastic strain, eep, by an increment, deep, that is proportional to de, 

deep — fde. 

The coefficient of proportionality, ip, may, in general, be a function of strain, e, stress, a, 
and the elastoplastic strain, eep. Only the dependence of </? on e is taken into account, 
which results in the kinematic equation 

^(i)=(^(e(t))|(t), eep(0)=0. (8) 

The function <^(e) vanishes at e = (the elastoplastic strain equals zero at very small 
strains), monotonically increases with strain, and reaches some limiting value b e (0, 1) at 
rather large strains (which corresponds to a steady regime of plastic flow) . To reduce the 
number of adjustable parameters in the constitutive equations, an exponential dependence 
is adopted, 

ip{e) =61 — exp(— ae) , (9) 

where the positive coefficients a and b are found by matching observations. 

Equations (8) and (9) differ from conventional flow rules in elastoplasticity, where 
the elastoplastic strain is assumed to be proportional to the stress, a. Similarities and 
differences between our approach and traditional ones will be discussed in Section 8. 
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5 Sliding of junctions in active MRs 



Sliding of junctions in active meso-domains with respect to their positions in a stress-free 
medium is treated as a rate-dependent process and is associated with the viscoelastic 
response of a semi-crystalhne polymer. Sliding of junctions in active MRs reflects (i) 
sliding of tie chains along lamellae, (ii) slip of chains in amorphous regions with respect 
to entanglements, and (iii) rearrangement of junctions driven by slip of lamellar blocks. 

The non-affine deformation of a network is modelled as a mechanically activated pro- 
cess induced by straining of active meso-domains. The strain in a meso-region, e, is 
defined as the difference between the macro-strain, e, and the elastoplastic strain, eep, 
caused by sliding of meso-domains with respect to each other, 

e(t) =e(i) -eep(i). 

Accepting the first-order kinetics for sliding of junctions, 

e{t) - eye{t,uj)], 
and using Eq. (1), we arrive at the constitutive equation 

eve(0,w) = 0, (10) 

where €ye{t,Lj) is the strain driven by sliding of junctions in an active MR with potential 
energy cu. 



■(i,u;)=r(a;) 



'-{t,uj) = roexp(-a;) e{t) - eep(t) - eye{t,(^) 



6 Strain energy density 

The elastic strain, ee, is calculated as the difference between the macro-strain, e, and the 
strains, eep and eye, induced by sliding of meso-domains with respect to each other and 
by sliding of junctions in active MRs with respect to their reference positions. 

For a passive meso-region (where sliding of junctions is prevented), the elastic strain 
is given by 

€e{t) = €{t) - eep(i). 
For an active meso-domain with potential energy cu, the elastic strain reads 

ee{t,uj) = e(t) - eep(t) - eve(t,a;). 

A strand is modelled as a linear elastic medium with the strain energy 

w = ^l^el, (11) 

where is a constant rigidity. 

Multiplying the mechanical energy of a strand, Eq. (11), by the number of strands per 
unit mass, summing the strain energies for strands in active and passive meso-domains. 
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and neglecting the energy of inter-chain interaction, we find the mechanical energy per 
unit mass of a semi-crystalline polymer, 



W{t) = -^^ 



r N,{t, u) [e{t) - eep(i) - e,,{t, uj)fduj + X^{t) (e{t) - e,^{t) 





(12) 



To develop a stress-strain relation, an explicit expression is necessary for the derivative 
of W with respect to time. Differentiation of Eq. (12) implies that 



dW 



it) 



fx 



X 



r N,{t,u)(^e{t) - e,p{t) - e,,{t,ij))du; + Xp{t)(e{t) - e,p{t) 

.J 



where 



"2^ 



too QN / \ 2 



+^it){e{t)-e^,{t)y 



roc , , Qf^ 

Y2{t) = /i N^{t,u;)[e{t)-e^p{t)-ey,{t,u;))-^{t,u;)duj. 

It follows from Eqs. (2), (3), (8) and (13) that 

dW 



dt 



it) = II 



X{e{t) - eep(i)) - P N^{t,uj)e^^{t,uj)duj 



1 - <f{eit))]f^{t) 



-Y,{t)-Y,{t). 
Equations (2) and (3) imply that 



ON, 



-{t, u})du}. 



dt dt Jo dt 

This equahty together with Eq. (14) results in 

Combining Eqs. (10) and (15), we find that 

Y2{t) = r N^{t,uj)T{tu)[e{t) - eep(t) - eve(t, c^))'o?u;. 

7 Constitutive equation 

For isothermal uniaxial deformation, the Clausius-Duhem inequality reads 

dW de 

— + a—>0. 

dt dt - 



duj. 



(13) 



(14) 
(15) 



(16) 



(17) 



(18) 



(19) 
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Substituting expression (16) into Eq. (19) and using Eqs. (4) and (5), we obtain 



a{t) - e(i - cfieit))) (e(t) - eep(t) - /^(t) p{t, uj)e,e{t, uj)duj) 
+Y,{t) + Y2{t)>0, 



de , , 



(20) 



where E = fiX. 

It follows from Eqs. (17) and (18) that for an active loading program (when e(t), e^dt), 
(t) and Nn{t,uj) increase with time), the functions Yi{t) and Y2{t) are non-negative. 
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This means that the dissipation inequality (20) is satisfied, provided that the expression 
in the square brackets vanishes, which results in the stress-strain relation 



a{t) = e(i - <p{e{t))) e{t) - e^p{t) - K{t) j%{t,u)e^^{t,u)du 



(21) 



Given functions K,(t), Q(t) and S(t), constitutive equations (6), (8) to (10) and (21) de- 
scribe the time-dependent response of a semi-crystalline polymer at uniaxial deformation. 
These relations are determined by 4 adjustable parameters: an analog of Young's modu- 
lus E in Eq. (21), dimensionless constants a and h in Eq. (9), and the attempt rate for 
sliding of junctions Fq in Eq. (10). The pre-factor Fq can be excluded from the governing 
equations by introducing a "shifted" potential energy uj = uj — ujq with uoq = IuFq. Be- 
cause this transformation does not change the structure of the stress-strain relations, we 
set Fq = 1 s without loss of generahty. 

For "rapid" deformations, when the effect of sliding of junctions in active MRs on the 
mechanical response of a specimen is negligible, the constitutive equations are substan- 
tially simplified. Neglecting the integral term in Eq. (21) and using Eq. (9), we find 
that 



E 



1 — b(l — exp(— ae 



where the elastoplastic strain eep obeys Eqs. (8) and (9), 



de, 



ep 



de 



b(l - exp(-ae)), eep(O) = 0. 



(22) 



(23) 



Equations (22) and (23) are determined by 3 material constants, E, a and b, to be found 
by matching a stress-strain curve for a tensile test with a constant strain rate. 

To study "slow" deformation processes, when sliding of junctions in meso-domains 
is to be taken into account, an additional hypothesis should be introduced to describe 
the effect of deformation history on the quantities k, Q and E. Two approaches are 
conventionally used to predict the effect of mechanical factors on these parameters, the 
so-called models with strain and stress clocks. For a survey of these concepts, the reader 
is referred to Drozdov (1998) and the bibliography therein. 

The theory of a stress-induced internal time is traditionally employed to fit observa- 
tions in creep tests with the loading program 



a{t) 



0, 



a 



t < 0, 
t > 0, 



(24) 
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where is a given stress. In terms of our model, this concept means that the quantities 
Q and E depend on the current stress a. Combining Eqs. (6), (8), (9), (21) and (24), 
we arrive at the formulas 



/•c 

'«(a°) / 
Jo 



2E2((70) 



+ 6^ 



(ie, 



cp 



dt 



it) 



l-6(l-exp(-ar")) 
^P^l-6(l-exp(-ae(i)))' 

6(l-exp(-ae(t)))^(t), eep(O) 



-ep- 



(25) 
(26) 



In these equations, the initial instant t = corresponds to the beginning of the creep test, 
the quantities e° and e[!p are found by integration of Eqs. (22) and (23) in the interval 
from o" = to 0" = o"°, the coefficient is determined by condition (7), and the function 
eyeit,uj) obeys Eq. (10). 

It is worth noting that Eq. (26) can be integrated exphcitly to express the elasto- 
plastic strain, eep, by means of the macro-strain e. However, we do not dwell on this 
tr ansf ormat ion . 

Given a stress, cr°, Eqs. (10), (25) and (26) are determined by 3 experimental con- 
stants, «;(cr°), n(cr°) and E((7°) to be found by matching observations in a creep test. 

According to the concept of a strain-induced material clock, the parameters k, Q, and E 
are functions of the current strain e. This approach is conventionally used to approximate 
experimental data in a relaxation test with 



e{t) 



0, 



t < 0, 
t > 0, 



(27) 



where e° is a given strain. It follows from Eqs. (8), (9) and (27) that the elastoplastic 
strain, eep, is time-independent. The quantity €ep{t) — e^p is determined by integration of 
Eq. (23) from zero to e°. Combining Eqs. (6), (21) and (27), we find that 



a{t) 



1-6 1 



exp(— ae°)j |f 

roo 

Jo 



-cp 



2E2(eO) 



(28) 



where the instant t = corresponds to the beginning of a relaxation test, the coefficient 
is determined by condition (7), and the function eye{t,ui) is governed by Eq. (10) with 

ro = i, 



dt 



expi^— a; J 



— — eve(t, tjj) 



eve (0,0;) 



0. 



Introducing the notation 



eve(^,o;) 



eve(t,a;) 



-ep 



we present the latter equation as follows: 
de 



dt 



{t,oj) = exp(-a;)[l - eye{t,oj)], eve(0,a;) = 0. 



(29) 



(30) 
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Substitution of expression (29) into Eq. (28) results in 



a{t) 



a 



.0 




2E2(eO) 




(31) 



where 



(7°(e°) ^E[l-b(l- exp(-ae°))] (e° - e°p). 



Given a strain, e'', Eqs. (30) and (31) are characterized by 4 experimental constants, 
<7^{e^), K{e^), fl{e^) and S(e°) to be determined by fitting observations in a relaxation 
test. 

Our aim now is to find adjustable parameters in the constitutive equations by matching 
experimental data depicted in Figures 1, 2 and 4, and to assess the apphcability of the 
concept of internal time with stress- and strain-induced material clocks. 

8 Fitting of observations 

We begin with the approximation of the stress-strain curve depicted in Figure 5. The re- 
striction on strains (cmax = 0.06) may explained by two reasons: (i) at higher elongations, 
the assumption that the strain energy of a strand is a quadratic function of strain, see 
Eq. (11), becomes questionable, and (n) according to Figure 3, at e = 0.06 the primary 
creep is transformed into the secondary creep (developed plastic flow), which is beyond 
the scope of the present study. 

Under uniaxial tension with the cross-head speed 5 mm/min, the strain emax = 0.06 
is reached within 69 s. According to Figure 2, changes in strain induced by sliding 
of junctions during this period are insignificant at stresses up to 20 MPa, whereas the 
duration of stretching at higher stresses does not exceeds 30 s, which causes rather small 
growth of strains. Based on these observations, we treat the deformation process as rapid 
and apply Eqs. (22) and (23) to match observations. To find the constants E, a and b, we 
fix the intervals [0, Omax] and [0, 6max], where the "best-fit" parameters a and b are assumed 
to be located, and divide these intervals into J subintervals by the points = iAa and 
bj = jAb (i, J = 1, .... J) with Aa = a^aax/ J and Ab = b^aax/ J- For any pair, {oj, bj}, we 
integrate Eqs. (22) and (23) numerically (with the step Ae = 5.0 ■ 10^^) by the Runge- 
Kutta method. The elastic modulus E = Eo{i,j) is found by the least-squares method 
from the condition of minimum of the function 



where the sum is calculated over all experimental points, ei, depicted in Figure 5, cxexp is 
the engineering stress measured in the tensile test, and (Jnum is given by Eq. (22). The 
"best-fit" parameters a and b minimize K on the set {ai,bj — 1, . . . , J)}. After 

determining their values, a, and bj, this procedure is repeated twice for the new intervals 
[ai_i,ai+i] and [bj-i,bj+i] to ensure an acceptable accuracy of fitting. The "best-fit" 
parameters read E = 2.12 GPa, a = 38.10 and b = 0.64. This value of E slightly exceeds 
Young's modulus {E — 1.50 GPa) provided by the supplier for a virgin material, which 
may be explained by changes in the microstructure of spherulites at annealing. 



= l][f^exp(e;) 
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To estimate the elastoplastic strain, egp, and the difference between the present model 
and the conventional flow rule 

(32) 

with a constant coefficient k, we present results of numerical simulation in Figure 6. This 
figure shows that the elastoplastic strain, egp, is neghgible at strains e < 0.025, and it 
increases (practically linearly) with macro-strain at higher elongation ratios. The ratio 

r(e) = ^ 

linearly grows with strain at e < 0.025, and it becomes practically constant at e > 0.03. 
The latter implies that in the region, where the influence of elastoplastic strains on the 
response of iPP cannot be disregarded, Eqs. (22) and (23) are rather close to the flow 
rule (32). An advantage of Eqs. (22) and (23) is that (i) they have a transparent physical 
meaning, (ii) no additional constants (analogs of the yield stress, (Xy, or yield strain, Cy) 
are to be introduced in the stress-strain relations, and (iii) these equations can correctly 
predict the time-dependent response of iPP in creep tests (see later), where Eq. (32) 
appears to be an oversimplifled relationship. 

Wc proceed with fitting observations in creep tests depicted in Figure 2. For any 
stress, (7°, the quantities e° and e^^ arc found by integration of Eqs. (22) and (23) with the 
material constants found in the approximation of the stress-strain curve plotted in Figure 
5. The quantities K(cr°), r2(cr°) and S(cr°) are determined by the following algorithm. We 
fix the intervals [0, Kmax], [0, f^max] and [0, Smax], where the "best-fit" parameters k, Q and 
E are assumed to be located, and divide these intervals into J subintervals by the points 
Ki = iAK, Qj = jAQ and = A;AS {i, j, k = 1, . . . , J) with Ak, = /tmax/ J, AQ = fimax/ J 
and AE = Emax/-/- For any pair, {Qj, Efc}, the constant = p^{j, k) is determined by 
Eq. (7), where the integral is evaluated by Simpson's method with 200 points and the 
step Ao; = 0.15. For any triple, {^j, Q^, E^}, Eqs. (10), (25) and (26) are integrated 
numerically (with the time step At = 0.1) by the Runge-Kutta method. The "best-fit" 
parameters k, Vt and E are found from the condition of minimum of the function 



K(i,j,k) = Xl[eexp(iO - e„um(iO 



where the sum is calculated over all experimental points, ti, presented in Figure 2, eexp is 
the strain measured in the creep test, and Cnum is given by Eq. (25). After determining 
the "best-fit" values, Ki, Qj and E^, this procedure is repeated for the new intervals 
[ftj-i, Kj+i], [flj_i,flj+i] and [Efc_i,Efc+i] to provide an acceptable accuracy of fitting. 
Figure 2 demonstrates fair agreement between the experimental data and the results of 
numerical simulation. 

The adjustable parameters Q, E and k arc plotted versus the engineering stress, a, in 
Figures 7 to 9. The experimental data are approximated by the linear functions 

Q = + ^iC, E = Eq -|- Ei(7, K — Ko + KiO", (33) 

where the coefficients Q^, E^ and {m — 0,1) are found by the least-squares technique. 
Figures 7 and 8 show that the quantities ft and E increase with stress, a, reach their 
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maxima in the interval between cr = 25 and cr = 30 MPa (which means, in the vicinity 
of the yield strain Cyi), and dramatically decrease at higher stresses. According to Figure 
9, the concentration of active meso-regions, k, increases with strain up to cr = 25 MPa, 
whereas at higher strains, the slope of the straightline (33) noticeably grows. 

Finally, we approximate the experimental data in the relaxation tests presented in 
Figure 4. To fit observations, we re- write Eq. (31) as 



a{t) = C,{e') +p'C2{e') J^^ e^t^oo) exp[-^^^^^^^]dw (34) 
with 

Ci = (7°, C2 = -(7°K. 

For any strain, e°, we fix the intervals [0, Jlmax] and [0, Smax], where the "best-fit" param- 
eters Q and S arc located, and divide these intervals into J subintervals by the points 
Qi = iAQ and = j AS {i, j = 1, . . . , J) with AQ = fimax/ J and AS = Smax/ J- For any 
pair, {fli, Sj}, the constant = p^{i,j) is determined from Eq. (7), where the integral is 
evaluated by Simpson's method with 200 points and the step Ao; = 0.15. Equations (30) 
are integrated numerically (with the time step At — 0.1) by the Runge-Kutta method. 
The constants Ci and C2 in Eq. (34) are determined by the least-squares algorithm to 
minimize the function ^ 

= X][o-exp(tz) - O-num(^0] > 
tl 

where the sum is calculated over all experimental points, ti, presented in Figure 4, Uexp is 
the stress measured in the relaxation test, and o"num is given by Eq. (34). The "best-fit" 
parameters Q and S are found from the condition of minimum of the function K{i,j). 
After determining the "best-fit" values, Qi and S^, this procedure is repeated for the 
new intervals and [Sj_i,Sj+i] to guarantee good accuracy of fitting. Figure 

4 demonstrates an acceptable agreement between the experimental data and the results 
of numerical analysis. 

The quantities Q, S and k = —C2/C1 arc plotted versus strain e in Figures 10 to 12. 
The experimental data are approximated by the linear functions 

Q = Qo + ^1^, S = So -|- Sie, K — kq + Ki€, (35) 

where the coefficients Qm, and k,„ (m = 0, 1) are found by the least-squares tech- 
nique. Figures 10 to 12 demonstrate that the adjustable parameters fl, S and k, found by 
matching observations in relaxation tests monotonically increase with strain e. 

There is no doubt that Eqs. (30) and (31) can be applied to fit experimental data in the 
relaxation test with the smallest strain = 0.05. Some questions arise, however, whether 
these equations (based on the quadratic approximation (11) of the strain energy), may 
be used to match observations in relaxation tests at higher strains. Fortunately, it can 
be shown that Eqs. (30) and (31) remain valid for an arbitrary (not necessary quadratic 
in strains) mechanical energy per strand w. To avoid complicated algebra employed in 
the derivation of constitutive equations at finite strains, appropriate transformations are 
omitted. 
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9 Discussion 



Our objective now is to compare adjustable parameters in the constitutive equations found 
by fitting observations in the creep and relaxation tests. 

First, we plot the quantities Q, E and k determined in the approximation of the 
relaxation curve at together with the data obtained by matching the creep curves 
(Figures 7 to 9). For this purpose, we replace the strain in the relaxation test, e?, by 
the corresponding stress at the beginning of the relaxation process, = Ci(e5). Figures 
7 and 8 demonstrate good agreement between the parameters Q and E determined by 
fitting the experimental data in the creep and relaxation tests. Figure 9 reveals that the 
concentration of active meso-domains, k, found in the relaxation text follows Eq. (33) 
with the coefficients found in the approximation of the creep curves below the yield point, 
Syi, see curve 1. These conclusions confirm changes in the distribution function p in 
the interval [eyi,ey2] revealed by matching the experimental data in the creep tests, but 
question appropriate findings for the concentration of active MRs k. 

To shed some light on this controversy, we re-plot the experimental constants fl, S 
and K found by fitting observations in the creep tests below the yield strain eyi together 
with those determined in the approximation of the relaxation curves (Figures 10 to 12). 
The adjustable parameters Q, E and k found by matching the creep curves are depicted 
as functions of strain at the beginning of creep tests. Figures 10 and 11 reveal that the 
quantities fl and S linearly increase with strain up to the yield strain eyi, where they suffer 
a pronounced drop, and proceed to grow with strain afterwards. Figure 12 demonstrates 
that the concentration of active meso-domains, k, linearly increases with strain in the 
entire region of strains under consideration, and the data found in the creep tests are in 
fair agreement with those determined by matching the relaxation curves. 

Based on these observations and adopting the dual lamellar population model (Verma 
et al., 1996), we propose the following scenario for the effect of mechanical factors on the 
distribution of active MRs: 

1. Below the yield strain Cyi, the average activation energy, fl, and the standard devi- 
ation of activation energies, E, for sliding of junctions in active MRs monotonically 
increase with strain e. This growth is attributed to fragmentation of thin (sub- 
sidiary) lamellae formed during injection-molding of specimens and developed at 
annealing of iPP. The increase in the average activation energy, fl, is associated 
with 

• breakage of subsidiary and transverse lamellae into small pieces that serve as 
extra physical cross-links in amorphous regions, 

• mechanically-induced activation of rigid amorphous fraction, where sliding of 
junctions was prevented by surrounding lamellae in stress-free specimens. 

The increase in the standard deviation of potential energies, E, reflects heterogeneity 
of the fragmentation process, which grows with strain because of the inhomogene- 
ity of breakage of subsidiary lamellae in active meso-regions with different sizes 
(potential energies). 
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2. When the strain, e, reaches the yield strain, eyi, subsidiary (thin) lamellae become 
totally disintegrated, which results in a decrease in the average activation energy 
Q (lamellar blocks that served as physical cross-links disappear) and a pronounced 
decrease in the standard deviation of potential energies S (meso-domains become 
more homogeneous). Figures 7 and 8 reveal that the interval of strains, where this 
"homogenization" of the micro-structure occurs, is relatively narrow (less than 2 
%), which imphes that the slopes of curves 2 in these figures are rather high. 

3. Straining of a specimen above the yield strain, Cyi, causes fragmentation of dominant 
(thick) lamellae, which results in a noticeable increase in Q and S. This growth is 
driven by the same mechanism as an increase in Q and E in the sub-yield region 
of deformation (breakage of lamellae into small pieces that serve as extra physical 
cross-links in the amorphous phase). This conclusion may be confirmed by the 
similarity of slopes of curves 1 and 2 for the standard deviation of potential energies 
of MRs, S, depicted in Figure 11. The slopes of the curves 1 and 2 for the average 
potential energy for shding of junctions, Q, in Figure 10 differ substantially. This 
discrepancy is attributed to the fact that the growth of the average potential energy 
of active meso-domains in the sub-yield region of deformation is governed by two 
morphological transformations (release of rigid amorphous fraction and formation 
of extra junctions in active MRs), whereas only the latter process takes place in the 
post-yield region. 

4. The concentration of active meso-domains, k, monotonically increases under stretch- 
ing. A noticeable increase in k. with stress reported in Figure 9 is an artifact caused 
by the assumption that the concentration of active MRs remains constant during a 
creep test (which, in turn, is based on a concept of stress-induced internal time). 
Despite good agreement between the experimental data and the results of numer- 
ical simulation revealed in Figure 2, it seems more adequate to presume that the 
adjustable parameters fl, S and n depend on the current strain (not stress). This 
hypothesis does not affect curves 1 to 4 in Figure 2 (because of relatively small 
increases in strain during the creep tests), but it results in substantial changes in 
creep curves 5 and 6 (corresponding to the yield region). Assuming k to be a func- 
tion of strain, one can treat the experimental points on curve 2 in Figure 9 as some 
"average" (over the creep curves) values of the concentration of active MRs in the 
post-yield region of deformation. In this case, the values of k found by matching 
creep curves 5 and 6 in Figure 2 become quite comparable with the data depicted in 
Figure 12, where these values roughly correspond to the vicinity of the yield point 
ey2. 

10 Concluding remarks 

A model has been developed for the elastoplastic and viscoplastic responses of semi- 
crystalline polymers at isothermal uniaxial deformation with small strains. To derive 
constitutive equations, a complicated micro-structure of a polymer is replaced by an equiv- 
alent permanent network of macromolecules bridged by junctions (physical cross-links. 
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entanglements and lamellae). The network is thought of as an ensemble of meso- regions 
with various potential energies for sliding of junctions with respect to their positions in a 
stress-free medium. 

The elastoplastic (rate-independent) behavior of a semi- crystalline polymer is attribut- 
ed to sliding of meso-domains with respect to each other. The kinetics of sliding is 
governed by Eqs. (8) and (9) with 2 adjustable parameters: a and b. The viscoplastic 
(rate-dependent) response is associated with sliding of junctions in active MRs. The 
sliding process is assumed to be thermally-agitated, and its rate is described by the 
Eyring equation (1). Equation (10) for the rate of viscoplastic strain is based on the first 
order kinetics of sliding, and it does not contain experimental constants (the attempt rate 
To is set to be 1 s). 

The distribution of active meso-domains is determined by the random energy model, 
Eq. (6), with three adjustable parameters: the average activation energy for sliding 
of nodes. Q, the standard deviation of potential energies of active MRs, E, and the 
concentration of active meso-domains, k, that are affected by mechanical factors. With 
reference to the concept of mechanically-induced material clocks, two hypotheses are 
analyzed: (i) the quantities Q, E and k are stress-dependent, and (ii) these parameters 
are governed by the current strain. 

Stress-strain relations arc derived by using the laws of thermodynamics. The consti- 
tutive equations are developed under the assumptions that (i) a strand can be modelled 
as a linear elastic medium with the strain energy (11), and (ii) the energy of inter-chain 
interaction is negligible. These equations contain 6 material constants: an elastic modulus 
E, and the quantities a, b, Q, E and k. 

To find these parameters, a tensile test with a constant strain rate, a scries of 8 creep 
tests, and a series of 4 relaxation tests have been performed on injection- molded isotactic 
polypropylene at room temperature. Material constants are determined by fitting exper- 
imental data. Figures 2, 4 and 5 demonstrate fair agreement between the observations 
and the results of numerical simulation. 

The following conclusions are drawn: 

1. Stretching of iPP in the sub-yield region of deformation increases the average poten- 
tial energy for sliding of junctions and the standard deviation of potential energies 
of active MRs. This increase is attributed to fragmentation of subsidiary lamellae 
and activation of the rigid amorphous fraction. 

2. In the vicinity of the yield point, the parameters Q and E fall down noticeably, 
which is associated with disappearance of physical cross-links formed by blocks of 
disintegrated thin lamellae. 

3. In the post-yield region of deformation, the quantities Q and E grow with strain, 
which is explained by fragmentation of dominant lamellae. 

4. The concentration of active meso-regions, K, monotonically increases with strain 
both in the sub-yield and post-yield domains. 
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5. The hypothesis about a strain-induced material clock appears to be more adequate 
to morphological transformations in isotactic polypropylene under loading than the 
assumption about the dependence of material parameters on the current stress. 

This study focused on modelling the nonlinear elastoplastic and viscoplastic behavior of 
iPP in creep and relaxation tests at relatively small strains. Some important questions 
remained, however, beyond the scope of the present work. In particular, transitions 
from the primary creep to the secondary and tertiary creeps have not been analyzed. 
No attention has been paid to the effect of time and temperature of annealing on the 
time-dependent response of isotactic polypropylene. The applicability of the constitutive 
equations to the description of the mechanical behavior of other semi-crystaUine polymers 
has not been examined. These issues will be the subject of a subsequent publication. 
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Figure 1: The engineering stress a MPa (unfilled circles) and the true stress Ut MPa 
(filled circles) versus elongation ratio A in a tensile test. Symbols: experimental data 

Figure 2: The strain e versus time t s in a tensile creep test with an engineering stress 
(T° MPa. Circles: experimental data. Curve 1: a\ = 10.0; curve 2: = 15.0; curve 3: 
= 20.0; curve 4: = 25.0; curve 5: cxg = 30.0; curve 6: = 30.38. Sohd fines: 
results of numerical simulation 

Figure 3: The Hencky strain en versus time t s in a tensile creep test with an engineering 
stress (7° MPa. Symbols: experimental data. Unfilled circles: 0-4 = 25.00; filled circles: 
aO = 30.00; asterisks: = 30.38; diamonds: = 30.94; triangles: = 32.80. The lines 

AA' and BB' indicate the strains corresponding to transitions from the primary creep to 
the secondary creep and from the secondary creep to the ternary creep, respectively 

Figure 4: The engineering stress a MPa versus time t s in a tensile relaxation test with 
a longitudinal strain e°. Circles: experimental data. Solid lines: results of numerical 
simulation. Curve 1: e? = 0.05; curve 2: = 0.10; curve 3: £3 = 0.15; curve 4: = 0.20 

Figure 5: The engineering stress a MPa versus strain e in a tensile test. Circles: exper- 
imental data. Solid line: results of numerical simulation 

Figure 6: The elastoplastic strain eep (curve 1) and the ratio r of the rate of elastoplastic 
strain to the engineering stress (curve 2) versus strain e in a tensile test. Solid lines: 
results of numerical simulation 

Figure 7: The average potential energy of meso-regions Q versus engineering stress a 
MPa. Symbols: treatment of observations. Unfilled circles: creep tests; filled circle: 
relaxation test. Solid lines: approximation of the experimental data by Eq. (33). Curve 
1: rio = 5.62, fii = 0.17; curve 2: = 25.66, = -0.60 

Figure 8: The standard deviation of potential energies of meso-regions S versus engi- 
neering stress cr MPa. Symbols: treatment of observations. Unfilled circles: creep tests; 
filled circle: relaxation test. Solid lines: approximation of the experimental data by Eq. 
(33). Curve 1: Eq = 2.52, Ei = 0.09; curve 2: Eq = 10.28, Ei = -0.26 

Figure 9: The concentration of active meso-regions k versus engineering stress a MPa. 
Symbols: treatment of observations. Unfilled circles: creep tests; filled circle: relaxation 
test. Solid lines: approximation of the experimental data by Eq. (33). Curve 1: Kq — 
0.0900, Ki = 0.0080; curve 2: kq = -1.1995, = 0.0595 

Figure 10: The average potential energy of meso-regions fl versus longitudinal strain e. 
Symbols: treatment of observations. Unfilled circles: creep tests; filled circles: relaxation 
tests. Solid fines: approximation of the experimental data by Eq. (35). Curve 1: = 
6.66, Qi = 130.66; curve 2: = 3.95, Qi = 20.20 

Figure 11: The standard deviation of potential energies of meso-regions E versus longi- 
tudinal strain e. Symbols: treatment of observations. Unfilled circles: creep tests; filled 
circles: relaxation tests. Solid lines: approximation of the experimental data by Eq. (35). 



22 



Curve 1: Sq = 3.53, Ei = 25.88; curve 2: Eq = 1.00, Si = 23.80 

Figure 12: The concentration of active meso-regions k, versus longitudinal strain e. 
Symbols: treatment of observations. Unfilled circles: creep tests; filled circles: relaxation 
tests. Solid line: approximation of the experimental data by Eq. (35) with kq — 0.16, 
Ki = 3.95 
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